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Abstract 


A  tmified  approach  to  the  theory  of  correlations  in  a  plasm 
is  presented,  based  on  the  3BKGY  hierarchy.  The  theory  is  applied 
to  a  one- component  plasma  with  the  Coulomb  interaction  modified 
to  include  effects  of  the  background.  Closed  integro-dif ferential 
equations  in  space  and  time  are  obtained  for  the  two-particle 
correlation  function  in  both  the  strong  and  weak  coupling  limits. 
In  the  weak- coup ling  domain,  \  <<  1,  the  time- independent  analysis 
returns  the  well-known  linearized  Debye-Huckle  result,  where  y  is 
the  plasma  parameter.  In  the  strong-coupling  domain  with  \  i  1, 
the  resulting  two-particle  'total'  correlation  function  exhibits 
decaying  oscillatory  behavior  for  particle  separation  of  the  order 
of  the  effective  interparticle  range.  .. 
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I .  Introduction 

Various  efforts  have  undertaken  to  describe  one- component  plasmas 

in  two  extreme  limits.  These  are  the  weakly  (y  <<  1)  and  strongly 

(y  i  1)  coupled  domains.  In  the  former  limit  many  such  efforts 

12  3 

have  led  to  the  linearized  Debye-Hiickle  result, using  kinetic  , 

statistical  mechanical^'^’^,  and  numerical^ ’ ® ^  formulations.  ! 

In  the  strong -coup ling  domain,  studies  have  centered  around  such 

problems  as  plasma  turbulence^^ ’ and  electron  correlations  at  metallic 
•  •  ''2 

densities  .  Each  approach  has  applied  its  own  approximation  of  the 
correlation  functions  to  the  BBKGY  hierarchy  in  order  to  obtain  the 
dielectric  response  functioiT  .  The  function  thus  acquired  is  then  em¬ 
ployed  in  conjunction  with  the  fluctuation-dissipation  theorem  to 
obtain  a  self-consistent  solution  for  the  static  form  factor  or,  equiv¬ 
alently,  the  pair  correlation  function"‘^. 

A  more  recent  study  and  review  of  calculations  concerning  the 

total  internal  energy  of  a  one-component  plasma  in  the  strong-coupling 

15 

limit  is  presented  by  Gould,  et  al.  Comprehensive  expositions  on 
the  state  of  the  art  of  this  subject  have  been  given  by  Kalman  and 
Carini;^^  Baus  and  Hansen,^  and  Ichimaru.^^ 

The  relevance  of  strongly-coupled  plasmas  to  natural  phenom.ena 
is  illustrated  in  Fig.  1,  We  note  that,  in  particular,  the  strong¬ 
coupling  limit  plays  an  important  role  in  the  description  of  x-ray 
plasmas,  laser  fusion  devices,  and  in  the  interiors  of  certain  super- 
dense  stars. 

In  the  numerical  work  of  Brush,  Sahlin  and  Teller',  the  pair 
correlation  function  is  evaluated  over  a  large  range  of  values  of 
the  plasma  paramter,  They  find  that  with  increasing  v  the 


f 


J 


correlation  function  passes  from  a  Debye-Hiickle  form  to  a  decaying, 
oscillatory  form.  This  nimerical  study  further  demonstrates  that 
the  effect  of  the  background  charge  is  to  alter  the  effective  inter- 
oarticle  interaction  away  from  the  bare  Coulomb  force. 

In  the  present  study  a  unified  formulation  for  the  theory  of 
correlations  in  a  one- component  plasma  is  introduced  which  is  valid 
for  weakly  and  strongly-coupled  plasmas.  The  theory  is  based  on  the 
33KGY  hierarchy  and  a  sequential  ordering  of  the  correlation  functions 
This  ordering  is  specified  in  terms  of  the  plasma  parameter  v.  An 
^-expansion  imbedded  within  the  ordering  scheme  permits  an  interactive 
technique  of  solution  and  renders  the  analysis  self-consistent. 

To  better  incorporate  the  role  of  the  neutralizing  background, 

the  effective  interparticle  force  is  expanded  about  the  bare  Coulomb 

interaction.  This  inclusion  is  motivated  by  results  of  previous 
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numerical  studies.  ’  In  both  the  weak  and  strongly  coupled 
domains,  closed  space-time  integro-dif ferential  equations  are 
obtained  for  the  two-particle  correlation  functions . 

In  the  weak-coupling  domain  a  time- independent  analysis  returns 

1_Q 

the  well-lmown  linearized  Debye-Hiickle  correlation  function. 

In  the  strongly-coupled  domain  a  second-order  differential 
equation  is  obtained  for  the  correlation  function.  The  solution  to 
this  equation  exhibits  a  decaying  oscillatory  structure  for  particle 
separations  in  excess  of  the  order  of  the  effective  two-particle 
interaction  range.  These  purely  analytic  results  are  in  very  good 
agreement  with  previous  numerical  work.' 

The  essential  components  of  this  analysis  are  related  as  shown 
in  the  flow  chart  given  in  Fig.  2. 


Analysis 


A.  Basic  Formulation 

We  consider  an  aggregate  of  charges  Ze  in  a  uniform  neutralizinj 
rH 

background.  The  s  equation  of  the  BBKGY  hierarchy,  hereafter  call^ 
BYg ,  is  given  by  (see,  for  example,  Liboff^^) 

r  (1) 

18 

Parameters  are  nondimensionalized  as  follows  ,  where  barred  variable 
are  dimensional. 


X  =  r^x,  D  =  o  p 
~  o~  -  o— 


G.  .  =  —  G.  ,  .  6 
-'ll  ~ij  ’  o 


-2  2 
Z  e 


In  addition, 

Pq  “  m  C  ,  m  C  =  kgT 

The  parameters  and  in  (1)  which  emerge  as  a  consecuence 
of  nondimensionalization  are. 


.  -  =  4Tnr^ 

O  O  S  V  o 

o 


Here  and  r^  are,  respectively,  the  strength  and  range  of  the  inter¬ 
action  potential.  For  the  case  that  r  =  the  Debye  distance,  the 

O  Cl 

coefficients  a  and  y  are  tnritten  without  subscriots , 


L 
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The  nondimens ional  distribution  is  then  given  by 

Fg  =  (mC)^®V®?g  (4) 

where  V  is  the  volume  of  the  system  and  is  the  s-oarticle  joint 
probability  distribution.  The  function  has  the  normalization, 

|Fg  dld2...ds  =  1  (5) 


where  phase-volume  elements  are  likewise  dimensionless.  The  effective 

til  t'l  “ 

force  on  the  i'^  ion  due  to  the  j  ion  is  . 

As  noted  in  the  introduction,  it  has  been  found  in  Monte  Carlo 
calculations^  that  the  net  effect  of  the  background  charge  is  to 
alter  the  effective  interparticle  force  away  from  a  simple  Coulomb 
form.  .‘Accordingly  we  write 


o  X.  . 


I 

n=l 


n,  2 
^n 


n-i 


(6) 


(6a) 


The  coefficient  u  is  a  parameter  of  smallness.  The  coefficients 
will  be  addressed  later  in  the  paper.  Note  that  the  leading  cor¬ 
rection  term  in  the  series  (6)  has  the  effect  of  altering  the  magnitude 
but  not  the  form  of  the  Coulomb  interaction. 

The  correlation  functions  are  defined  through  the  Mayer  expansions 


F2(1,2)  =  F^(l)  F^(2)  ^  C2(1,2) 

F3(1,2,3)  =  Ft^CD  F3_(2)  F^(3)  (7 

^  F3_(1)  C2(2,3)  +  F^(2)  C2(3.1)  ^  F^(3)  Cod.l) 


03(1,2,3) 


One  also  writes 


02(1,2)  =  Fj^Cl)  7^<i2)  h(l,2)  (8) 
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where  h(l,2)  is  the  so-called  'total  correlation  function' 

In  the  weak -coup ling  (or,  'correlationless')  domain,  one  assumes 

^1^1  C2  >>  .  (9a) 

By  extension,  in  the  strong-coupling  domain  we  write 

Fj_F^  <<  <<  C3  ...  (9b) 

These  sequential  inequalities  may  both  be  incorporated  into  the 
single  set  of  expansions, 

Fi  =  y°[F^^°^+£F^^^V  ...] 

C2  =  yCC^°^  +  s  +  ...]  (10) 

C3  =  +  e  +  ...] 


The  £ -factors  follow  standard  perturbation  procedure  and  permit  an 
iterative  technique  of  solution.  In  (10)  we  have  reintroduced  the 
plasma  parameter, 


which  in  this  instance  serves  as  a  bookkeeping  index.  With  v  <<  L,  I 

(10)  is  seen  to  return  the  weak-coupling  sequence  (9a) ,  whereas  wirh 
V  >  1,  (10)  yields  the  strong-coupling  sequence  (9b).  Thus  the  expansion 
(10)  serves  as  a  unified  description  of  plasma  conditions.  The  Debye 
distance  and  plasma  frequency  which  enter  (11)  are  given  by 


4T(eZ)' 


k  T 
n 


2.,2 

:ne  Z 


m 


.  2  2 
‘d  '"p 


7 


B .  The  Strongly  Coupled  Domain^ 

In  this  limit  we  take  -x^  ->  1,  and  the  relation  y  *  gives 
v.  Thus  BY 2^  and  BY2  become,  in  accord  with  expansion  (7), 


'1^1  =  -  +  c^] 


(K2  +  aQB2)CF2_(l)F2_(2)  +  C2]  = 

2  . 

-  J-  IjCF,  (1)F,  (2)F,  (3)  +  :  Ft  (i)C-,(j  ,k) 

4tt  z  i  i  i  P(ijk)  ^ 


(15a) 


+03(1.2.3)]  (15b) 

Here  the  symbol  P(i,j ,k)  denoces  summation  over  permutations  of  ijk. 

Our  chief  aim  at  this  point  is  to  obtain  a  closed  kinetic 
equation  for  €2^°^  (1,2).  To  this  end  a  study  was  made  of  the 

c-dependence  of  the  parameters  v  and  c ,  with  '^1.  It  was  found 
2  -1/2  2 

that  Y  e  and  uy  1  are  the  simplest  forms  which  give  closure, 
(See  Appendix  A.)  Applying  the  expansions  (10)  and  (14)  to 
(15)  with  these  £ -dependencies  of  v  and  „  gives  the  following 
leading  perturbative  equations  relevant  to  determining 


-  'll  = 


(16a) 


-  h  Ii^°^5'i(l)Fi(2)  =  0 


(16b) 


(:<2  I^(-)  :  F^(i)c.^°^:  ,k)  ( 

2  0  2  .  1 

Here  we  have  written  F^^  for  We  note  in  passing  that  (17) 

is  a  closed  space-time  kinetic  equation  for  Co ^  ^(1,2).  Con¬ 
straints  (i5a,b)  are  consistent  with  an  assumption  of  spatial 
homogeneity,  so  that  we  may  write 


7.1 


The  right-hand  side  of  (17)  is  simplified  by  virtue  of  (16a) 
and  the  vanishing  of  any  space  integral  over  an  isotropic  vector 
field.  Consequently  two  terms  survive  given  by 

..2b  2 

RHS(17)  =  -  -2^  {F^(2)  il;-  +F^(1)  ^2)  •  £2 


[Go  to  p.  8] 
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where 


=  Jd3  (x,3) 


(19  a) 


(19b) 


With  reference  to  the  results  of  Appendix  B,  we  may  write 


(dropping  the  superscript  on  h) 

T  (1.0)  ,  ^  -2  r 

J  '  ^  3  -  4ttx  X  / 

"  y  C 


g  dc  3  h(p) 


,  (1.0)  _  ,  (1.0) 

~2  ■  ■  tl 


(20a) 


(20b) 


Here  and  henceforth  we  set  x  h  x^2 • 

With  (17)  and  (19)  in  (16)  we  obtain 

(Ko  +  3^B2'^°^F,^(P3_)F3_(p3)h(x) 


(21) 

X  2 

do  o‘'h(p) 


In  the  space-dependent  part  of  K2  we  use  the  transformation 


Xi  +  X2  =  X  x^  -  X2  =  X 


in  which  case 


K2h(x)  =  X  •  (p^  -  P2)  -r^  h(x)  (22) 

where  we  have  neglected  the  time  derivative  in  K2.  With  the  repulsive 
bare  Coulomb  interaction,  it  follows  from  (5)  and  (14)  that 

*2'°^  ■  ?  •  ■  '4] 

Employing  the  preceding  two  expressions  in  (21)  gives 
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X 


K  2  2 

bi  Y 


W  ■  4l]  ''l<Pl)^l<P2>'^ 

4:  ■  4]  ^i'pi^-i^p2>  •  ? 

(24) 


Assvnning  finally  chat  is  isotropic  in  momentum  space,  we  may 
write 


g 

3p 


^l^P)  “  P  ^^1^P> 

and  (24)  becomes 

3 


(25) 


X  •  p. 


-2^1 ^P2^  3 

^1^1  lx  ■’■  ^o^  ~p^^  ^l^Pl^ 


-b,  Y 


2  2  ^1^P2^  3 
'W~  5Pl 


-^(p^)x-2  f 

J  o 


dp  P  h(p) 


(26) 


^  r  "  -2  ^i^pi^ 

?  •  52  i;  *  “o’'  ^  ’'l(P2’ 


2  2  ^l^Pl^ 


-Vv‘^^PlCP2>--‘f  P‘«,) 

operating  on  this  equation  with  J" dV2  eliminates  the  RHS,  and  since 
/\ 

X  •  p^  is  generally  nonvanishing,  there  remains 

''l^Pl>  ^  ^  5P^  'l^Pl^ 


SP2 


-  b 


^  -±  w  X 


X 


Pi  '^Pl 

Integration  over  p^  gives 


h(x)  -  Y  Y“x~^h(x)  -i-  b 

where 


do  o^h(o)  =  0 


2  2„2  -2  f 

Jo 


do  o^h(o)  =  0 


(23) 


i 

1 
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^  -/■'Pi  ^  it 


and  may  be  assijmed  positive. 

Differentiating  (28)  with  respect  to  x  gives 
2a 

h"(x)  +  —2 —  h(x)  -  a^Y^x'^h'Cx) 


+  b,^Y^Y^  I  h(x)  - 


do  0  h(p) 


Using  (28)  to  solve  for  the  integral  term  in  (30)  yields  the 
differencial  equation 


h"  + 


[2y’" 


'7  J  h'  +  h  =  0 


where 


y  s  kx  k  =  b^Yy  <  = 


Ola) 


23 

Employing  the  transforaation 

h(y)  =  1  e"^^^u(y) 


a  =  </2 


in  (31)  gives 


u"  +  (1  -  ^)u  =  0  (33) 

7 

24 

This  equation  is  Schroedinger-like  in  form  and  we  may  conclude 
that  for  y“  >  a,  h(y)  is  oscillatory,  whereas  ror  y“  <  r,  h(y)  is 
e:cponencial.  The  explicit  form  of  the  solutions  in  these  domains  may 
be  obtained  in  wo  extremes : 


V  >>  a 


We  are  left  with 
u”  +  u  *  0 

which  has  the  solution 


u  =  .4  cos  V  +  3  sin  V 


With  (32)  we  obtain 


>  e  '  r  T 

h  (y)  »  — - -  [A  sin  y  +  B  cos  yj 


In  this  case  (33)  reduces  to 
2 

u”  -  2^  u  -  0 

y 

24 

which  has  the  solution 


(35) 


(36) 


u(y)  =  y(A  -f  B  e"'''^^) 

With  (32) ,  this  gives 

h^(y)  -  B  e”"^y  (37) 

where  we  have  set  A  =  0  because  (31)  does  not  allow  a  constant 

solution.  As  previously  noted,  the  oscillatory  behavior  (35)  of 

h(y)  in  the  domain  y  >>  a,  is  in  very  good  agreement  with 

7  8 

previous  numerical  studies  '  . 

We  now  estimate  the  value  of  the  leading 

coefficient  b^  in  (6) .  This  is  accomplished  through  comparison  of 

the  wave  number  K  (defined  below)  of  h^(y)  with  the  numerical  value 

7  8 

obtained  in  Monte  Carlo  studies  '  .  These  results  imply  that  in 
the  vicinity  of  y  =  50  (i.e.,  r  =  10) 

“  If  - 


3  =  1.5  .  r 

The  wavenumber  K  is  given  by 
y  =•  Kic 


K  =• 


8.4 


(38) 


From  (11),  Furthermore,  assuming  a  Maxwellian  form  for 


/ 


Fi(pi)  gives  Y  “  1.  With  these  values 
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K  -  h^/X^ 


and  comparison  with  yields 

which  with  the  values  cited  above  gives  “  0.8.  With  this  result 
at  hand,  we  may  estimate  the  value  of  interparticle  separation  at 
which  oscillation  of  h(x)  ensues.  Such  oscillation  was  found 
from  (33)  to  occur  for 

>  </2 

or,  equivalently 


X 

r_ 


2b^ 


Inserting  our  previous  finding 


b-,  Y  =  lY 


.2/3  ^  2/3 


gives 


X  >  r^/Y 


1/3 


Thus  oscillation  of  h(x)  may  be  expected  for  interparticle 
displacement  greater  than  or  equal  to  the  range  of  the  two-particle 
interaction.  The  initiation  of  such  oscillation  may  be  interpreted 
as  the  onset  of  phase  transition  in  the  plasma  ’  Furthermore 
the  wavenumber  K  of  these  oscillations  (33) ,  is  seen  to  grow  with 
Y.  This  behavior  may  be  associated  with  the  fluid-solid  phase 
transition^  at  large  y,  found  in  numerical  studies.’’ 
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C.  The  Weak-Coupling  Case 

This  sicuation  has  been  studied  by  many  individuals^”^  and 
as  we  have  noted  previoxisly  it  gives  the  linearized  Debye-Hiickle 
result.  In  the  kinetic  domain  this  limit  yields  the  Vlasov  equa¬ 
tion,  or  more  generally,  the  Bales cu- Lenar d  Equation. 

An  essential  element  of  these  studies  is  that  the  correlation 
functions  are  assumed  to  be  perturbativelv  small.  Accordingly  we 

set  Y  -  £  in  (10).  Insofar  as  r^  =  X,  in  this  limit,  it  follows 

o  Q  ’ 

2 

that  n  =  ~j./y  ®  1,  which  fixes  :i  ~  .  We  note  that  this  limit  is 

O  ■? 

equivalent  to  the  classical  Rosenbluth-Ro stoker  ansatz.  With 
these  constraints  and  u  ~  (6)  gives  the  following  leading 

equations  for  €2^*^^  as  derived  from  (13) 


Vi  '  ^  ^1  ^1  ^1 


(39a) 


(39b) 


iCiFi^^^d)  =  -  ^  I],^°^CC2^°^(1,2)+F2*^°^(1)F2‘'^^(2)+F^^'^^(1)F2^°^(2)  j 

K2[C2^°^  (1,2)+Fi‘^°^  (l)Fi^^^  (2)+Fi^*\l)Fi‘^°\2)  1  +  iB2^°^F2^°\l)  F^^®^  (2)  = 


(39c) 


“  F  <°)(i)c/°^(jk)+  : 

^  ,P(ijk)  ^  “  P(ijk) 


p(o)(i)F(o)(j)F(l)(k) 


We  must  solve  (39)  and  (40)  simultaneously  for  C.,^°^(l,2). 
First  note  that  (39a)  may  be  rewritten  in  two  equivalent  ways: 

F^^^^(2)K^(1)F^'^°^(1)  =  -  ^  Fj_‘^"^(2)l2‘^°^(l)F2*^°^(l)F^'^°^(3) 
F^^^^(1)K2(2)F2'‘^°^2)  =  -  ~  F;^'^^l)l2^°^(2)F,^^°^(2)F/'^\3) 


13.1 


Similarly  (39b)  gives 

Fj^''°^(2)Kj^(1)f/^^  (1)  -  -  (2)1^^°^  (1)  (1,3)  + 

F^^°^(1)F^^^^(3)  +  Fj^^^^(l)Fj^<^°^(3)] 

F^^°^(1)K^(2)F^^^^(2)  -  -  ^  F^^°^(1)I^^‘^^(2)[C2^°^(2,3)  + 

Fj_^°^  (2)F^^^^  (3)  +  F^^^^  (2)F^^°^  (3)  ] 


Adding  these  four  equations  and  subtracting  the  result  from  (39c) 
yields 


K2C2^°^(1.2)+aB2*^°^F^<^°^Fj_^°^  =  -  ^  -^2 P  (i]k)  ^1 ^2 


Fj_^°^  (2)I^^°^(1)C2^°^  (1.3)-F^^°^  (1)I3_'^°^(2)C2^°^(2,3)  :■ 


Four  terms  remain  on  the  RHS  of  this  equation.  Of  these,  two  vanish 
by  integration  of  the  isotropic  vector  field  ,  leaving 


K.,C2^°\i,2)  +  = 

_2 


(40) 


■  i;  {Ii_‘'°^(l)Fj^^°^(l)C2^°^(2.3)  +  l2(2)Fj^*^°^(2)C2^°^  (1,3)  } 

Here  we  have  obtained  (as  in  the  strongly-coupled  analysis)  a 
closed  space- time-dependent  equation  of  notion  for  0.,^°^.  The 
RHS  of  this  equation  simplifies  to 

RHS(40)  =  -  ^  F^^°^  (1)  - 


-  (o) 
‘1 


(1) 


'Po 


F,  (2) 


(41) 


13.2 


where  Che  J  integrals  are  defined  in  (19)  ,  and  again  we  find  that 


T  (0.0) 


(42a) 


T  (0.0)  ^  ,  (0.0) 

ii  ■  .2 


(42b) 


[Go  to  p.  14] 
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Substituting  (18)  and  (42)  into  (40)  we  obtain 


K2FiFih(x)  +  aB2^°^Fj_Fj_ 


-  n 


I 


Fj^Fj^  •  X  X 


-2 


; 

/  do  p 

Jo 


■b(o)  (43) 


Making  use  of  (22)  and  (23)  in  the  LHS  of  (43)  gives 


3S  *  ^?1  "  Ba^^l^l^  +  ax  •  X 

2  r 

-  n  I 


-2 


"81  382 


^1^1 


^1^1 


-  -2  n 

X  X 

Jo 


L^Bi  ^B2j 

Furthermore,  with  (25)  the  preceding  equation  becomes 

a  ^i^F2^ 


do  p  h(p) 


(44) 


X  . 


Pi 


F,F,  X  h  +  F,  (p,) 

1  1  3X  ^  p^  jp^  l^^l^ 


2  ^1^P2^  -2  r 

+  n  ^  F3_(Pj^)x  ^  / 

n 


?l  5P;j_ 


dp  p  h(p) 


(45) 


S  •  P2 


3  a  ^l^Pl^  ' 

^1^1 


2  ^l^'Pl^  a  -2  2 

^  ^^^F^(P2)x  2  j  do  o2h(p) 


P2  JP2  -  -  -o 

Repeating  the  procedure  leading  to  (27)  ,  we  find  that  (45) 
yields 


Integration  over  p-,  gives 


^  h  -  4  -  p  V--’ 

x“ 


o 


X 


■h(p)  =  0 


("463 


(47) 


*> 


ifferentiating  (47)  leads  to 


h-  +  ^  +  2n  Vx'^ 


do  o^h(p)  -  n^Y^h  *  0 


Combining  the  last  two  equations  gives 


h"  + 


h  -  0 


where 


7  =  r,Yx  (49a) 

The  general  solution  to  (48)  may  be  written 

h(7)  -  + 

We  need  to  affirm  that  (49)  is  a  solution  to  the  parent  integral 
equation  (47) ,  which  in  terms  of  y  is 

d  1  r  ^  /'y  2 

^  h  -  -j  otY'^n  +  /  dy  y'^h(y)  =0  (51) 

y  L  Jo 

Inserting  (49)  into  (50)  leads  to  the  constraint 

Di  +  D2  =  -  aY^n  (52) 

Setting  D2  -  0  for  physical  reasons  leaves  the  solution 


h(y) 


aY-^n  -y 

s  •  a  v 


For  a  Maxwellian,  Y  =  1.  Furthermore,  in  the  present  limit  we 
recall  that  n  =  1,  so  that  y  =  x,  and  (52)  may  be  written 

h(x)  =  -  (54) 

35 

This  is  the  well-known  linearized  Debye-Huckle  result  ,  which 
follows  to  first-order  in  a  from  an  expansion  of  the  more  precise 
nonlinear  Debye-Hiickle  form 

h(x)  -  -  1  (55) 


where 
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(55a) 

Whereas  the  Debye-Huckle  form  (54)  agrees  with  the  linearized  result 
(53)  for  X  >>  1,  we  see  that  at  the  origin  h(x)  diverges,  but  h^  has 
the  correct  value  h^(o)  =  -  1. 


III.  Conclusions 

We  have  described  a  unified  formulation  of  correlations  in  a 
one-component  plasma,  based  on  the  plasma  parameter  y,  which  is  valid 
in  both  the  weakly  and  strongly-coupled  domains.  The  theory  is 
based  on  the  BBKGY  hierarchy  together  with  expansions  of  the  two- 
particle  interaction  and  correlation  functions. 

In  the  weak-coupling  domain,  the  analysis  returns  the  familiar 
linearized  Debye-Hiickle  correlation  function.  In  the  strong- coup  ling 
domain  a  second-order  differential  equation  for  the  total  correlation 
function  is  found  which  yields  an  exponentially  damped  solution  near 
the  origin.  For  particle  separation  on  the  order  of  the  range  of  the 
two-particle  interaction,  the  equation  gives  a  Bessel-like  oscillatory 
solution.  These  features  of  the  correlation  function  are  in  excel¬ 
lent  agreement  with  previous  numerical  studies. 

In  both  domains,  closed  space-time  dependent  integro- 
differential  equations  were  obtained  for  the  two-particle  correla¬ 
tion  functions.  We  anticipate  that  future  study  will  find  impor¬ 
tant  applications  of  these  results  both  in  plasma  and  condensed 
matter  physics. 


Anoendix  A 


Here  we  examine  the  decoupling  condition  on  y  employed  in 
Section  3.  With  expansions  (6)  and  (10)  in  BY^  (15b),  we  write 

°  t=o  ^  ^  n=0  ^  J 


(Al) 


2  r »  /  ^ 

1 _  r  ^  *  f  p  T?  7  J.  .. 

.  !  -  j-o  .r-ii-ir, 

^■'  ir=0  2  I  1  -  ^ 

L  J  ^ 


,n|^(n)  j_  .  2 


Z  _mp  (m)"] 


The  ?T  functions  are  expanded  as  in  (10).  We  seek  a  kinetic 
equation  for  consistent  with  aQ*0(l).  In  general,  such  an 

/V  ^  ,v 

equation  includes  K,  3  and  I  terms.  The  mere  requirement  of  a  K 
term  indicates  that  the  equation  sought  must  be  0(y) »  which  leads 
to  the  condition 

^  T*  3 

Y  =  -4“','  =  Y  (A2) 

The  leftmost  equality  gives  t  =  0,  with  u  1 .  This  specifies  inclu¬ 
sion  of  the-  32^^  term  in  our  equation.  The  rightmost  equality  gives 


The  simplest  choice  xvhich  includes  effects  of  the  background  is 


r  =  1  which  soecifies  inclusion  of  the  I 


(1) 


term  ana 


o 

V 


(A3) 


Finally,  decoupling  the  term  from  ail  other  terms 


■2  z 

[in  the  RHS  of  (Al)]  reauires  that 


0  n 


0  ,  n 

0  ,  m 


(A4) 

(A5) 
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These  relations  give,  respectively, 


q  i  0  ,  n  5 

1 

(A6) 

S-1  !Ii  ,  , 

YU  £  #  1 

s  ^  0  ,  mi 

0 

(A7) 

Using  (A3)  in  (A6)  yields 

^  ^n/2(q-l) 

q  i  0  ,  n  y 

1 

(A8) 

The  strong-coupling  condition  (y  ^  1)  requires 

,n/2(q-l)  ^  ^  ^ 

q  i  0  ,  n 

5  1 

(A9) 

This  allows  as  the  only  possibility,  q  =  0.  T 

herefore , 

(A8) 

requires 

-n/ 

Y  ?  £ 

0 

n  u  1 

(AlO) 

Siailarly,  with  (A3)  in  (A7)  it 

follows  chat 

.m/2s-3  , 

C.  >  1  , 

s  u  0  ,  mi 

0 

(All) 

This  leads  to  the  cases  s  =  0,l 

and  thus 

,,  ^  .-n/3 
{  r  t 

£*“  m  >.  0 

(A12) 

Combining  results  (AlO)  and  (Al 

2) ,  the  final 

conditions 

on 

!  are 

summarized  by 

-a/b 

V  =  £  a  u  1  ,  0 

u  4  ,  f  ^  m 

m 

’  2  ’ 

m 

3 

(Ail) 

where  n  5  Q . 
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Any  one  of  these  allowable  orders  of  y  leads  to  the  sought-after 
decoupled  equation  for  in  the  strong-coupling  liait.  In  this 

analysis  we  select  the  simplest  case:  a  =  l,  b  =  4,  corresponding 
to 


2 

Y 


,-1/2 


(A14) 


[Go  to  p.  20] 
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Aooendix  B 


Consider  the  integral  (19a) .  We  set 
£  -  ?23  ^ 

Xi3  =  3  +  X 

in  which  case  for  n  =  1,  m  =  0,  we  find 

T  (1.0)  _  3  r.-  ^(h) 

£].  ^  J-  15  +  X,  <B2 

In  spherical  coordinates  with 
2 

dc  =  0  sin  0  d0  d<p  do 

0  •  X  =  xp  cos  9  (33 

!o  +  xl  =  (p^  +  2x*o  x^)^ 

and  X  held  fixed  along  the  z-axis  throughout  the  integration  we 


obtain 


j^d.O)  ^  2^  J  dp  p^h(p)  J  dS  sin  0  (p^  +  2xpcos  3  +  x^)"^ 


Performing  the  3-integration  gives 


-  -  2’ 

-  0 

or.  equivalently, 


dp  p^h(p) 


-  -  2^  ik  /J  d-- 


[o^  +  2xpcosG  -t-  x-J  “ 


?  9  if 

p^  -  2x0  +  x“)^  -  (o"  *  2x0  +  X 


2’  A  r  f  h(c)[;.-  - 

-  ^  o 


X  '  -  !  0  4-  X  ' 


A  #  W.)  .  fji.  ch(.) 
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Figure  Captions 
Figure  1 

The  plasma  parameter  y  n  and  T  (with  2=1)  showing 
weakly  and  strongly  coupled  (shaded)  domains. 

Figure  2 

Flow  chart  explaining  the  analysis. 

Figure  3 

ct^-y  plot  showing  weakly  coupled  (WCP)  and  strongly 
coupled  (SCP)  plasma  domains. 
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V^alues  of  the  plasma  parameter  y  as  a  function  of  electron  number  densit\' 
and  temperature. 
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